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ABSTRACT 

Among the icy satellites of Saturn, lapetus shows a striking dichotomy between its 
leading and trailing hemispheres, the former being significantly darker than the lat- 
ter. Thanks to the VIMS imaging spectrometer on-board Cassini, it is now possible to 
investigate the spectral features of the satellites in Saturn system within a wider spec- 
tral range and with an enhanced accuracy than with previously available data. In this 
work, we present an application of the G-mode method to the high resolution, visible 
and near infrared data of Phoebe, lapetus and Hyperion collected by Cassini/VIMS, 
in order to search for compositional correlations. We also present the results of a dy- 
namical study on the efficiency of lapetus in capturing dust grains travelling inward 
in Saturn system with the aim of evaluating the viability of Poynting- Robertson drag 
as the physical mechanism transferring the dark material to the satellite. The results 
of spectroscopic classification are used jointly with the ones of the dynamical study 
to describe a plausible physical scenario for the origin of lapetus' dichotomy. Our 
work shows that mass transfer from the outer Saturnian system is an efficient mecha- 
nism, particularly for the range of sizes hypothesised for the particles composing the 
newly discovered outer ring around Saturn. Both spectral and dynamical data indicate 
Phoebe as the main source of the dark material. However, due to considerations on 
the coUisional history of the Saturnian irregular satellites and to the differences in the 
spectral features of Phoebe and the dark material on lapetus in the visible and ultra- 
violet range, we suggest a multi-source scenario where now extinct prograde satellites 
and the disruptive impacts that generated the putative coUisional families played a 
significant role in supplying the original amount of dark material. 

Key words: planets and satelUtes: general; techniques: spectroscopic; methods: nu- 
merical; planets and satellites: individual: lapetus; planets and satellites: individual: 
Phoebe; planets and satellites: individual: Hyperion 



1 INTRODUCTION 

Saturn's third largest satellite, lapetus, shows one of the 
most striking dichotomies in the solar system. Observations 
of the satellite, starting from a couple of years after its 
discovery by G.D. Cassini in 1671 and culminating in two 
Voyager flybys and several Earth-based investigations (e.g., 
Cruikshank et al. (1983); Squyres & Sagan (1983); Squyres 
et al. (1984)), had shown that the trailing hemisphere has 
an albedo, composition, and morphology typical of the other 
icy Saturnian satellites, being highly reflective and spectrally 
consistent with water ice, while the surface of the leading 
hemisphere is composed of a very low surface reflectance 
(geometric albedo 2-6%), reddish material that may include 
organics and carbon. 
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The origin and nature of the low-albedo hemisphere of lape- 
tus is still debated, representing one of the most intriguing 
problems in planetary science. The two main theories of the 
origin of the low-albedo hemisphere are that it was created 
by an endogenic geologic process such as flooding by mag- 
mas (Smith et al. 1981), or that it resulted from accretion 
of exogenous particles (Soter 1974; Matthews 1992). The 
primary evidence for the endogenic model is the existence of 
in-fllled craters (Smith et al. 1981), while tlie main support- 
ing observation for the exogenous model is that the dark ma- 
terial is centred precisely on the apex of motion (Matthews 
1992). 

Disk-integrated spectra collected by the VIMS instrument 
on-board Cassini during several serendipitous periods in 
2004 July and October (at distances between 2-3.1 x 10® 
km) enabled a first analysis of lapetus' surface in the spec- 
tral range from 1 to 5 /im that led to the identification of 
CO2, which is most prominent on the dark side (Buratti et 
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al. 2005). These authors also found that a good fit to the 
low-albedo average spectrum is obtained by a small amount 
of ice, an even smaller amount of Fe2 03 to account for the 
ferric absorption band at ~1.0 /im, 36% Triton tholin, and 
substantial amounts of the HCN polymer ( "poly-HCN" ) . 
On 31 December 2004, the Cassini spacecraft had a distant 
flyby with lapetus (> 124,000 km, yielding a maximum spa- 
tial resolution on the satellite's surface of ~62 km), that 
represented the first good chance to acquire spatially re- 
solved spectra of the satellite with better signal-to-rioisc ra- 
tio (SNR) in the range beyond 3 /im, though in this CEise the 
relatively large phase angle (54° to 116°) allowed to explore 
only a fraction of the daysidc. where most of the framed 
surface was covered by dark material. On these data, Cruik- 
shank et al. (2008) reported the identification of a broad 
absorption band centred at 3.29 /jm, and concentrated in a 
region comprising about 15% of the low-albedo surface area, 
that they interpreted as the C-H stretching mode vibra- 
tion in polycyclic aromatic hydrocarbon (PAH) molecules. 
Cruikshank et al. (2008) also found, in association with the 
aromatic band, two weaker bands at 3.42 and 3.52 /um, at- 
tributed to the asymmetric stretching modes of the -CH2- 
group in aliphatic hydrocarbons. Compounds containing the 
nitrile (-C=N) functional group were initially suspected as 
the cause of the 2.42-2.44 ^m band, clearly correlated to the 
low-albedo material on Phoebe, lapetus and Dione (Cruik- 
shank et al. 1991; Clark et al. 2005; Brown et al. 2006; 
Clark et al. 2007), and possibly some of the indistinct spec- 
tral structure in the region 4.5-4.8 /Ltm (Buratti et al. 2005); 
although recently Clark et al. (2009) revised this interpre- 
tation and suggested that molecular hydrogen trapped in 
dark material would be a better candidate for this feature. 
From the analysis of far ultraviolet spectra returned by the 
Cassini Ultraviolet Imaging Spectrograph (UVIS) , Hendrix & 
Hansen (2008) found that water ice amounts increase within 
the dark material away from the apex (at 90 W° longitude, 
the centre of the dark leading hemisphere), consistently with 
thermal segregation of water ice; yet the fact that water ice 
is present also at the lowest, darkest and warmest latitudes, 
where it is not expected to be stable, may be a sign of on- 
going or recent emplacement of the dark material from an 
exogenous source (ibid). 

Images returned by the Imaging Science Subsystem on- 
board Cassini during this flyby suggested mass wasting of 
ballistically deposited material (Porco et al. 2005), although 
the limited spatial resolution did not allow to safely identify 
any small bright crater lying over the dark terrain, which 
is a strong argument in favour of deposited material. More- 
over, using scatterometry data returned by the Cassini Ti- 
tan RADAR Mapper (named RADAR for concision in the 
following), Ostro et al. (2006) pointed out that lapetus' 2.2- 
cm (13.78 GHz) radar albedo is dramatically higher on the 
optically bright trailing side than the optically dark leading 
side, whereas 12.6-cm (2.38 GHz) results reported by Black 
et al. (2004) with the Arecibo Observatory's radar system 
show hardly any hemispheric asymmetry and give a mean 
radar reflectivity several times lower than the reflectivity 
measured at 2.2 cm, which can be explained if the leading 
side's optically dark contaminant is present to depths of at 
least one to several decimetres. 

On September 10, 2007, the Cassini spacecraft performed its 
first and only targeted lapetus flyby, at minimum altitude 



of 1620 km. Approach occurred over the unlit leading hemi- 
sphere, departure over the illuminated trailing one. From 
the ISS data, a strong indication that dark material is lying 
over bright terrain finally came from the numerous small 
bright craters detected in the high-resolution images from 
the flyby. The existence of these small craters (up to some 
dozens of meters in diameter) indicates that the dark blanket 
is very thin, probably no more than a few meters, possibly 
only decimetres. On the other hand, the dark material is 
expected to be thicker than only millimetres because of the 
evidence for mass wasting and because of the ability of the 
Cassini 2.2 cm RADAR to distinguish between the dark and 
the bright terrain. 

The confirmation of the exogenous model emphasised the 
need to identify a reliable source of the dark material. In 
past times, for different reasons, a number of possible sources 
were proposed: Phoebe (Soter 1974), Titan (Owen et al. 
2001), Hyperion (Matthews 1992), or possibly other small 
dark irregular satellites (Buratti et al. 2002; Buratti, Hicks 
& Davies 2005), although other works essentially focused 
on Hyperion and Phoebe as the main candidates. The mech- 
anisms invoked to justify the transfer of material from the 
source to lapetus were different for these two bodies: for 
Hyperion the proposed mechanism was coUisional excava- 
tion and/or break-up (Matthews 1992; Marchi et al. 2002; 
Dobrovolskis & Lissauer 2004), while Poynting- Robertson 
{PR in the following) drag was advocated for Phoebe (Soter 
1974; Burns et al. 1996). PR drag was also invoked by Bu- 
ratti, Hicks & Davies (2005) in their suggestion that the 
source of the dark material could be linked to the retrograde 
irregular satellites as a whole instead than just to Phoebe. 
Dust ejected from retrograde satellites and migrating inward 
from the outer Saturnian system would impact lapetus on its 
leading hemisphere due to its counter-revolving motion, thus 
naturally explaining the leading-trailing asymmetry (Soter 
1974). Intuitively, such mechanism should be characterised 
by a high transfer efficiency since lapetus' orbital period is 
shorter than the migration timescale by orders of magni- 
tude, yet to present the only quantitative evaluation of the 
transfer efficiency is that of Burns et al. (1996), who esti- 
mated it being about 70% for 10 /im sized dust grains. On 
the contrary, several evaluations of the transfer efficiency 
from Hyperion to lapetus have been performed (Matthews 
1992; Marchi et al. 2002; Dobrovolskis & Lissauer 2004), 
the resulting value varying over three orders of magnitude 
(from 0.1% to about 20-40%) depending on the characteris- 
tics of the ejection process. 

In addition to the limited information on the transfer effi- 
ciency, the major issue of the Phoebe-based model is that the 
visible spectrum of lapetus' leading side is similar to D-type 
asteroids, a primitive, very low albedo group of bodies ex- 
hibiting a typically reddish spectrum in the visual region, i.e. 
increasing in reflectance with increasing wavelength, while 
Phoebe's spectrum is essentially flat and grey in the visual 
region, thus similar to C-type or F-type asteroids (see for 
example Tholen & Zelhier (1983); Buratti et al. (2002); 
Grav & Bauer (2007)); moreover its albedo is higher than 
carbon, particularly in the brighter areas detected by Voy- 
ager 2 and Cassini. 

In this paper, we have explored VIMS data returned by 
the Cassini spacecraft after its close encounter with lape- 
tus, combining them with VIMS data of Phoebe and Hy- 
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perion - acquired with comparable spatial resolution and 
favourable illumination conditions - that were extensively 
analysed in previous works. We classify this dataset, with 
an automatic statistical method, separately for the visual 
and IR portions, to identify homogeneous types and look 
for correlations among the spectra of lapetus, Phoebe and 
Hyperion. To support the results of the spectral comparison, 
we reviewed and extended the dust migration scenario based 
on Poynting Robertson drag in light of the recent results of 
both theoretical works and the analysis of Galileo, Cassini 
and Spitzer data. We estimated its mass transfer efficiency 
by means of a statistical approach derived from the one de- 
signed by Kcsslcr (1981) to evaluate the impact probabilities 
of Jovian irregular satellites. Finally, we discuss the amount 
of ejecta needed to supply the dark material present on lape- 
tus and compare it to the one needed in the Hyperion-based 
scenario and the available observational constrains. 



2 THE G-MODE METHOD 

The G-mode method was originally developed by A.I. 
Gavrishin and A. Coradini (see Gavrishin, Coradini & 
Fulchignoni (1980); Gavrishin et al. (1992); Coradini, 
Pulchignoni & Gavrishin (1976); Coradini et al. (1977)) 
to classify lunar samples on the basis of the major oxides 
composition. The good results obtained warranted its appli- 
cation to several different data sets (sec, for example, Cora- 
dini, Fulchignoni & Gavrishin (1976); Carusi & Massaro 
(1978); Gavrishin, Coradini & Fulchignoni (1980); Bianchi 
et al. (1980); GiovanncUi et al. (1981); Coradini, Giovan- 
nelli & Polimene (1983); Barucci et al. (1987); Orosei et al. 
(2003)). In particular, the Imaging Spectrometer for Mars 
(ISM), flown on-board the Soviet Phobos mission, offered 
the first chance to apply the G-mode method to imaging 
spectroscopy data (Coradini et al. 1991; Erard, Ccrroni & 
Coradini 1991; Cerroni & Coradini 1995). More recently, 
the G-mode was applied also to Cassini/ VIMS data acquired 
during the close flyby of Phoebe (Tosi et al. 2005; Coradini 
et al. 2008) and to the combined data of Titan returned by 
Cassini/VIMS and Cassini/RADAR (Tosi et al. 2009). 
The G-mode differs from other broadly used statistical 
methods - such as the Principal Components Analysis 
(PCA) and the Q-rnode method - in some key characteristics 
(see Bianchi et al. (1980)): in summary, a linear dependence 
of the variables is not needed; instrumental errors can be 
taken into account; meaningless variables are discerned and 
removed; and finally different levels of classification can be 
performed. 

Basically, by lowering the confidence level of the test, set 
a priori by the user, the algorithm can perform a more re- 
fined classification, in order to look for further homogeneous 
types. In this case, the G-mode includes a test that allows 
to interrupt the classification when it becomes too detailed: 
when the statistical distance among types becomes smaller 
than the established confidence level, the algorithm can ei- 
ther stop or continue by merging different small types to- 
gether (this condition is reported in the output of the pro- 
gram) . 

In this specific case, our statistical universe is composed by 
N spectra sampled in M spectral channels. The G-mode al- 
lows the user to apply either the same error value to all the 



variables (this approach is preferred when all the variables 
shall drive the classification at the same level), or to assign 
a specific error value to each variable. In the case of multi- 
spectral and hyperspectral data, when using refiectances as 
variables, a logic approach is to use the instrumental noise, 
that for each spectral channel is given by the inverse of the 
SNR, as an error. Alternatively, should the depths of some 
diagnostic absorption bands be used as variables, then the 
average standard deviation of the band depth can be used 
as an error for each variable. In this work, we actually apply 
the first approach to the visible portion of the data, where 
no absorption band is expected to occur, while the second 
approach is preferred for the infrared data, where several 
spectral signatures are used to classify the data. 



3 THE VIMS INSTRUMENT 

The Visual and Infrared Mapping Spectrometer (VIMS) is an 
imaging spectrometer on-board the Cassini Orbiter space- 
craft. VIMS is actually made up of two spectrometers, 
VIMS-V (developed in Italy) and VIMS-IR (developed in 
the USA) . VIMS is the result of an international collabora- 
tion involving the space agencies of the United States, Italy, 
France and Germany as well as other academic and indus- 
trial partners. 

The two channels share a common electronics box and are 
co-aligned on a common optical pallet. The combined optical 
system generates 352 two-dimensional images (with maxi- 
mum nominal dimensions of 64x64, 0.5 mrad pixels), each 
one corresponding to a specific spectral channel. These im- 
ages are merged by the main electronics in order to produce 
"image cubes" representing the spectrum of the same field 
of view (FOV) in the range from 0.35 to 5.1 fim, sampled in 
352 bands. See Brown et al. (2004); Miller et al. (1996) for 
a complete description of the instrument. 
In this work, we separately explored both the visible and 
infrared portions of the selected cubes. 



4 THE DATASET 

Among the VIMS acquisitions, we selected 3 cubes, respec- 
tively from the Phoebe flyby and the to-date closest flybys 
with lapetus and Hyperion. Such data were selected in or- 
der to have: a spatial resolution as high as possible, a low 
phase angle (intended to select only cubes showing most of 
targets' dayside, and to limit the dependence on the spe- 
ciflc observational geometry), and a good SNR. To do this, 
we have selected cubes from the three satellites showing at 
the same time: 1) a spatial resolution ^ 3 km/pixel; 2) a 
phase angle < 43°; 3) an IR exposure time ^ 120 ms/pixel. 
Details about these data are given in Tables 2 and 3. The 
VIMS cube selected for Phoebe is the closest available from 
the flyby, acquired at 19:32 UT on 11 June 2004 at a range 
of 2096 km and phase angle 24.6°. Also for Hyperion we 
selected the closest cube, acquired at 02:15 on 26 Septem- 
ber 2005 at a range of 2989 km and a phase angle of 42.6°. 
Finally, the cube selected for lapetus was acquired at 15:01 
on 10 September 2007, at a range of 6091 km and a phase 
angle of 17.3°: although this cube was not acquired at clos- 
est approach but only 40 minutes later, it represents a good 
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Table 1. VIMS specifications summary. 







VIMS-V 




VIMS-IR 








Spectral coverage (/^m) 


0.35 - 1.05 




0.85 - 5.1 








Spectral channels (bands) 


96 




256 








Average spectral sampling (nm/channol) 


7.3 




16.6 








Total FOV (°) 


1.83x1.83 




1.83x1.83 








Total FOV (mrad) 


32x32 




32x32 








Nominal IFOV (mrad) 


0.50x0.50 




0.50x0.50 








Hi-res IFOV (mrad) 


0.167x0.167 




0.25x0.50 








Nominal image dimension (pixel) 1x1, 6x6, 12x12, 64x64 1x1, 6x6, 12x12, 64x64 






Detector type 


Si CCD (2D) 


InSb photodiodes (ID) 






Average instrumental SNR 


380 




100 






Table 2. Details of the VIMS cubes used in this work. 












Satellite 


Mission Dataset Filename 


Date 


UTC 


Relative 


Spacecraft 


Phase 




Sequence 




on-board 


velocity 


altitude 


angle 










(km/s) 


(km) 


n 


Phoebe 


SOI PHOEBE017 CM_1465674952 


11 Jun 2004 


19:32 


6.35 


2096 


24.576 


Hyperion 


S14 HYPERIONC007 CM_1506393701 


26 Sep 2005 


02:15 


5.64 


2989 


42.598 


lapetus 


S33 ORSHIRESOOl CM_1568129671 


10 Sep 2007 


15:01 


2.35 


6091 


17.298 
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Figure 1. RGB images (R=2.02 ^m, G=1.59 ^m, B=1.28 /im) of the VIMS cubes used in this work, shown in their IR portion. From 
left to right, the data refer to Phoebe (cube CM_1465674952), Hyperion (cube CM_1506393701) and lapetus (cube CM_1568129671), 
respectively. A horizontal scale bar, placed under each image, accounts for the linear dimension of the image. 



Table 3. Details of the VIMS cubes used in this work. 



Satellite 


Image 
dimension 

(PX) 


(msec/px) 


IR Spatial 
resolution 
(km/px) 


VIS tcxp 

(msec) 


VIS Spatial 
resolution 
(km/px) 


Phoebe 


30x18 


180 


0.524x1.048 


5630 


0.349x0.349 


Hyperion 


36x24 


320 


0.747x1.494 


7680 


0.498x0.498 


lapetus 


40x40 


120 


3.045x3.045 


4800 


1.015x1.015 
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Figure 2. RGB images (R,=702.9 nm, G=549.5 nm, B=439.2 nm) of the VIMS cubes used in this work, shown in their VIS portion. From 
left to right, the data refer to Phoebe (cube CM_1465674952), Hyperion (cube CM_1506393701) and lapetus (cube CM_1568129671), 
respectively. A horizontal scale bar, placed under each image, accounts for the linear dimension of the image. 



tradeoff among image dimensions, spatial resolution, expo- 
sure time and phase angle. Moreover, this cube is peculiar 
because it shows a region located exactly at the boundary 
between dark and bright terrains, where a relatively sharp 
division can be seen. In all cases, there are no sky back- 
ground pixels since the satellite's surface always fills the 
VIMS FOV. Figure 1 shows the RGB images of these cubes 
in their VIS portion. 

Due to the fact that the Phoebe and Hyperion cubes were 
acquired in high resolution IFOV mode both by the visible 
and infrared channels of VIMS, the IFOV of VIMS-IR being 
250x500 ^irad wide while the IFOV of VIMS-V is 167x167 
prad wide, in the case of these two satellites the IR and VIS 
images do not match. On the contrary, the lapetus cube 
was acquired in nominal IFOV mode (500x500 /xrad) by 
both channels, so that the IR image basically matches with 
the VIS image (actually a slight misalignment occurring be- 
tween the optical axes of the two VIS and IR telescopes 
prevents a perfect match). Figure 2 shows the RGB images 
of these cubes in their IR portion. All of the VIMS cubes 
used in this work were initially calibrated, in their infrared 
portion, by means of the RC15 VIMS-IR sensitivity func- 
tion and the flat-field cube released in 2005, namely the last 
official products available at the time this work was under- 
taken. The sensitivity function allows to convert the raw 
signal of each pixel inside the IR image into radiance (di- 
vided by the IR integration time and the fiat-field) and then 
into refiectance I /F, where / is the intensity of reflected light 
(uncorrected for the specific observational geometry and the 
thermal emission) and ttF is the plane-parallel flux of sun- 
light incident on the satellite (Thekaekara 1973), scaled for 
its heliocentric distance (for details about the VIMS cali- 
bration, see McCord et al. (2004)). It should be noted that 



the latest unofficial RC17 sensitivity function of VIMS-IR, 
derived in 2008 and currently under test, is different from 
the RC15 calibration used in Clark et al. (2005) by the use 
of additional standard stars and data from the VIMS so- 
lar port. It was found that the RG15 calibration, based on 
comparison to Galileo NIMS data with VIMS data from the 
Cassini Jupiter encounter, the Moon flyby, plus ground cal- 
ibration data contained a residual 3-micron absorption and 
spectral structure in the 2-2.5 micron region due to ringing 
in the VIMS order sorting fllters (Clark, personal commu- 
nication, 2009). For this reason, in order not to mislead our 
work, we have separately tested both these sensitivity func- 
tions on the VIMS infrared data. 

Moreover, all the spectra were "despiked": we removed 
single-pixel, single-spectral channel deviations, caused by 
systematic instrumental artifacts like order-sorting fflters as 
well as random artifacts like cosmic rays or high energy radi- 
ation striking the detectors (these events reveal themselves 
as spikes in the dark current stored in each raw cube) . On the 
visible portion of the VIMS data, also a "destriping" proce- 
dure was applied, aimed to remove residual offsets between 
radiometric levels of different CCD columns while observing 
a uniform scene. 

For the Phoebe and Hyperion cubes, we selected all the 
pixels inside the spectral image, so we have 540 samples 
from Phoebe and 864 samples from Hyperion. In the case 
of the lapetus cube, showing both dark and bright mate- 
rial, we deflne a "region of interest" including 639 samples 
from the dark material, covering 40% of the image. Hence 
it follows that, by summing the samples from Phoebe, la- 
petus and Hyperion, our global dataset is made up by 2043 
samples / spectra. 
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5 DISCUSSION 

5.1 Results in the visible range 

In the range measured by VIMS-V (0.35 - 1.05 ^m), no 
signatures can be safely identified in tlie spectra of lape- 
tus, Plioebe and Hyperion; hence a spectral classification 
is driven by the spectral continuum, namely the slope pa- 
rameter. First, all the 2043 calibrated spectra of lapetus, 
Phoebe and Hyperion were normalised to the value of 1.0 at 
the wavelength of 549.54 nm (VIMS channel 28), then the 
G-mode was applied with a 97.13% confidence level (1.90ct). 
The classification returned 5 homogeneous types, composed 
by 505, 744, 175, 538 and 71 samples respectively. 
A convenient way to represent the results is plotting the 
data in a multidimensional space with axes corresponding 
to the VIMS spectral channels nearest to the filters defined 
in the most common photometric systems: U (366.3 nm), B 
(446.5 nm), R (659.1 nm), / (805.2 nm) and Z (885.9 nm); 
the V filter being not considered, as it has been used as a 
reference wavelength for the normalisation of all the spec- 
tra. In this way, the three major types, i.e. type 1, 2 and 4 
(in blue, cyan and yellow, respectively), can be easily iden- 
tified in a 3D space whose axes correspond to the U, B and 
R photometric variables (see Fig. 3); for better clarity, the 
same data points can be also plotted in a 2D space with axes 
corresponding to pairs of the U, B, R, I and Z bands, so that 
the colour indexes for each sample can be evaluated as well: 
since we keep the information about the real object (satel- 
lite) corresponding to each sample, it can be verified that 
the main three types correspond to the three satellites (see 
Fig. 4) . Type 1 (blue) is made up by samples of the dark side 
of lapetus, that appear rather concentrated in all the photo- 
metric variables. Type 2 (cyan) is linked to Hyperion, show- 
ing the steeper slope of all three satellites (and thus being 
redder than the dark side of lapetus, as inferred by the high 
reflectance values in the R, I and Z bands) and with samples 
that show a larger spread with respect to the samples of la- 
petus. Type 4 (yellow) is clearly related to Phoebe, whose 
samples, despite of the small area of the satellite considered 
for this work, are more scattered in the multidimensional 
space with respect to the samples of lapetus and Hyperion 
(especially towards the shorter wavelengths) and are defi- 
nitely not reddened but, on the contrary, show subsets of 
pixels exhibiting a higher refiectance towards the blue and 
UV wavelengths, consistently with recent Earth-based ob- 
servations (Grav & Bauer 2007). Type 3 (green) includes a 
relatively large number of samples that show an intermedi- 
ate behaviour between the dark side of lapetus and Phoebe 
or Hyperion. Finally, type 5 (red) is made up by samples of 
Hyperion that behave differently from the majority of the 
samples of that satellite, showing higher reflectance values 
on an average in the R, I and Z bands and thus being redder 
than all the other types. 

This classiflcation points out the fact that the colour in- 
dexes of the dark side of lapetus, Phoebe and Hyperion, 
evaluated on high resolution data, are significantly different 
from each other. As expected, Hyperion is spectrally red at 
visible wavelengths, similar to D-type asteroids, with a sub- 
set of samples showing the steeper spectral slope of all the 
dataset; yet this reddening is, on an average, higher than the 
dark material of lapetus, while Phoebe represents a totally 
different case, with no reddening and even a minor number 



Phoebe 
lopetus 
Hyperion 






Type 


1 






2 




Type 


3 




Type 


4 
















Figure 3. Classification in the visible range. The samples are 
plotted in a 3D space, with axes corresponding to the normalised 
reflectances measured in the U, B, and R bands. Upper plot: 
distribution of the samples as a function of the three satellites. 
Lower plot: distribution of the homogeneous types returned by 
the G-mode analysis. 



of samples exhibiting a negative slope (i.e. the reflectance 
increases towards the shorter wavelengths). 

5.2 Results in the near infrared range 

In the range measured by VIMS-IR (0.85 - 5.1 fj,m), unlike 
the visible range case, several absorption features can be 
identified in the spectra of the three icy satellites. In the 
past, for classification purposes, we used the spectral range 
offered by VIMS-IR, both as a whole and truncated at 4.3 
^m in order to discard the contribution of thermal emission 
for lapetus, using simple refiectances I/F as variables (see 
Tosi et al. (2006)). In this work, we adopted a different ap- 
proach: we normalise all the spectra to VIMS channel 179, 
corresponding to the 2.23 wavelength, that is free from 
absorption features in all the three satellites, then we select 
a number of diagnostic absorption bands to be used as vari- 
ables. 

All the spectra from lapetus, Phoebe and Hyperion are dom- 
inated by the large 3 fxm absorption band due to the OH 
fundamental stretch in H2O ice and/or bound water, which 
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Figure 4. Normalised reflectances of the samples from Phoebe, Hyperion and the dark side of lapetus, projected on the UB, RB, RI, 
and IZ planes, respectively, from top to bottom. Left column: association with the satellites. Right column: homogeneous types found 
by the G-mode analysis; a black triangle marks the barycentre of the classes. 
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is therefore not relevant for the classification. Instead, we 
use the spectral signatures of H2O at 1.51 and 2.05 /im, 
combined to a number of spectral features diagnostic for 
volatiles. 

Interestingly, the us asymmetric stretch of CO2 is slightly 
shifted in position, moving - on an average - from 4.247 
Hm on Hyperion to 4.254 /im on lapetus and 4.260 ^m on 
Phoebe (Clark et al. 2005; Buratti et al. 2005; Cruikshank 
et al. 2007; Filacchione et al. 2009), i.e. a maximum dif- 
ference of 13 nm, which is anyway comprised into one spec- 
tral channel of VIMS-IR (whose spectral resolution at those 
wavelengths is ~22 nm). We sample this compound at 4.26 
/xm, the absorption being broad enough to prevent a signif- 
icant error to be committed in evaluating the band depth. 
The aromatic C-H stretch on Phoebe mostly occurs at 3.25 
/xm (Clark et al. 2005; Coradini et al. 2008) while on la- 
petus is reported at 3.29 fim (Cruikshank et al. 2008), i.e. 
shifted by ~40 nm or two spectral channels of VIMS-IR. 
We also include the CH aliphatic stretch at 3.41 fj.m and its 
overtone at 1.75 nm (Cruikshank et al. 2008); and the 2.44 
broad feature possibly being the overtone of a cyanide com- 
pound (Clark et al. 2005) or the fundamental of molecular 
hydrogen trapped in dark material (Clark et al. 2009). Ac- 
cording to the first interpretation, the fundamental of CN 
varies in position as a result of variations in its chemical 
environment; assuming the RC15 calibration, in the Phoebe 
spectrum a feature is seen in the 4.5-4.6 fim range (Clark 
et al. 2005) and we evaluate this bond at 4.53 fim. Finally, 
assuming the RC15 calibration we have also considered the 
4.42 /im signature, reported by Coradini et al. (2008) in 
homogeneous types identified on the surface of Phoebe and 
suggeted to be related to a nitrile compound like IIC3N or 
HNCO. It should be noted that by applying the RC17 cali- 
bration, the 4.53 pm and 4.42 pm features tend to disappear: 
such a behaviour can be indicative of calibration artifacts, 
so we discarded these two variables in our second test. 
The list of the spectral features used as variables in our clas- 
sification is summarised in Table 4. For each spectral signa- 
ture, we compute the band depth following the definition by 
Clark & Roush (1984), as: 



D = l 



Rb 

Rc 



(1) 



where Rb is the refiectance measured at the band centre 
and Rc is the reflectance of the spectral continuum at the 
band centre, reconstructed through a linear flt relying on 
the wings of the band. 

We also combine these band depths to the normalised re- 
fiectance of the spectral continuum sampled in two wave- 
lengths that are safely free from absorptions (at 1.82 /im 
and 3.55 /im), as a set of convenient variables to be used for 
the classification (see Table 4). 

By composing such a dataset and classifying it through 
the G-mode method with a 98.98% confidence level (2.05o-), 
we find 2 types made up of 849 and 1194 samples, respec- 
tively. By plotting the result in a 3D space with axes corre- 
sponding, for example, to the abundances of H2O ice (mea- 
sured at 1.51 /im), CO2 (measured at 4.26 /im) and C=N 
(measured at 4.53 /im), these types can be highlighted and 
the corresponding samples can be put in association with 
the satellite (see Fig. 5). 



Phoebe 
lapetus 
-Hyperion 






Figure 5. Classification in the near infrared range. The samples 
are plotted in a 3D space, with axes corresponding to the band 
depths of water ice (1.51 /im), CO2 (4.26 fim) and the CN funda- 
mental stretch (4.53 fim). Upper plot: distribution of the samples 
as a function of the three satellites. Lower plot: distribution of 
the homogeneous types returned by the G-mode analysis. 



For clarity, the same points can be also plotted in a 2D 
space with axes corresponding to the depths of some diag- 
nostic spectral features (aromatic/aliphatic CH, CO2, C=N, 
etc.) versus the 1.51 /im signature, related to the abundance 
of water ice. Keeping the information about the satellite con- 
nected to each sample, it can be easily verified that samples 
of type 1 (blue) entirely correspond to spectra of Hyperion, 
while samples of type 2 (red) are - with the exception of 
only one sample - spectra of Phoebe and of the dark side of 
lapetus (see Figs. 6 and 7). 

From the classification, it turns out that Hyperion and 
Phoebe are generally richer in water ice than the lapetus 
dark terrain, consistently with the results obtained in the 
far ultraviolet wavelengths by Cassini/UVIS (Hendrix & 
Hansen 2008), with a subset of the samples of lapetus and 
Phoebe exhibiting similar amounts of ice. The abundance of 
CO2 is more prominent in the lapetus dark material than 
on Phoebe, in agreement with the results by Buratti et al. 
(2005), with Hyperion representing an intermediate situa- 
tion and showing a stronger correlation with water ice. The 
abundances of hydrocarbons (traced by the CH aromatic 
and aliphatic stretches) are quite similar on the three bod- 
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Table 4. List of spectral features used to classify the VIMS data in the infrared range. 
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ies, but with a stronger correlation with water ice in the CEise 
of Hyperion. On the other hand, the dark material of lape- 
tus and Phoebe show similar amounts of nitrile compounds 
(assuming CN-bcaring molecules as traced by the 2.44, 4.42 
and 4.53 /im signatures), that however turn out again to be 
deeply bounded to water ice on Hyperion, whereas such a 
correlation is not found on the other two bodies. 
Regarding the spatial variability, for the considered dataset 
it should be noted that CO2 is more scattered in samples 
from lapetus and Hyperion than in samples of Phoebe. Ni- 
triles traced at 4.53 fj,m and 4.42 /im show a larger vari- 
ability on Hyperion, although the opposite is interestingly 
observed for the 2.44 /im feature, especially on Phoebe. The 
aromatic CH stretch is more variable on Hyperion, whereas 
the aliphatic CH signature shows a larger variability on la- 
petus. Finally, the spectral continuum sampled at 1.82 /xm 
(once the spectra have been normalised at 2.23 /im) turns 
out to be similar for Phoebe and lapetus, while Hyperion 
shows a higher I/F at that wavelength; the three satellites 
show anyway different I/Fs at 3.55 /um. 

After this analysis, we redid the work in the infrared 
range by applying the latest, unofficial sensitivity function 
derived for VIMS-IR (RC17, Roger Clark, personal commu- 
nication, 2008). As mentioned before, with this new calibra- 
tion the features at 4.42 and ~4.5 /xm disappear on most 
data, so they are no longer considered as variables for the 
classification, while we keep all the other features, so that 
the new dataset is made up by 2043 spectra in 9 variables. 
By classifying it through the G-mode method with the same 
confidence level used before (98.98% or 2.05cr), we find 2 
types made up of 844 and 1199 samples, respectively. Fig- 
ures 8 and 9 show the result of the classification, with sam- 
ples represented in a 2D space and axes again correspond- 
ing to the depths of the considered spectral features with 
respect to the 1.51 /im signature related to the abundance 
of water ice. It can be seen that the classification remains 
essentially unchanged with respect to the application of the 
RC15 sensitivity function, with the same variability of the 
key compounds on the satellites and a general strong corre- 
lation with water ice that still characterizes the samples of 
Hyperion in comparison to those of lapetus and Phoebe. 
In light of this analysis, a general evidence is that, regard- 
less of the specific sensitivity function used to calibrated 
the VIMS data, the situation in the near infrared range is 
different with respect to the visible range: on the basis of 
the depths of either nine or seven spectral features, related 
to the abundances of several compounds or spectrally ac- 



tive bonds on the surfaces of the satellites, and of the spec- 
tral continuum profile sampled in two different wavelengths, 
the existence of two classes defines a general higher degree 
of similarity between the dark hemisphere of lapetus and 
Phoebe; the G-mode method allows to come to this con- 
clusion through an automatic statistical approach applied 
to a large number of samples, thus increasing the degree of 
confidence. However, despite the small spatial scale of the 
considered data, this analysis points out a general stronger 
correlation with water ice for Hyperion, and a certain degree 
of variability of spectral features related to hydrocarbons, 
CO2 and other compounds on the three satellites. 



6 ORIGIN OF THE DUST: PHYSICAL AND 
DYNAMICAL SCENARIO 

As we mentioned in the introduction, Phoebe and Hyper- 
ion have historically been proposed as the main potential 
sources of the dark material coating the leading hemisphere 
of lapetus (Soter 1974; Matthews 1992). More recently, 
Buratti, Hicks & Davies (2005) suggested that the source 
of this material could be linked to more than one object, 
namely the retrograde Saturnian irregular satellites. To de- 
livery the dark material from Hyperion to lapetus-crossing 
orbits, dynamical models relied on high ejection velocities 
(generally v > 100 m s"S see Matthews (1992); Marchi et 
al. (2002); Dobrovolskis & Lissauer (2004)) and possibly 
on collimated clouds of cjecta (Matthews 1992; Marchi et 
al. 2002). PR drag, causing the dust grains to spiral inward 
toward Saturn, was instead invoked by those scenarios link- 
ing the dark material to Phoebe and the irregular satellites 
(Soter 1974; Buratti, Hicks & Davies 2005). 
In this section we will discuss the dust production and trans- 
fer mechanisms in the outer Saturnian system (i.e. the region 
populated by the irregular satellites) and compare them with 
the ones estimated for the Hyperion-based scenario. 



6.1 Irregular satellites: impacts and dust 
production 

The images of Phoebe taken by Cassini mission (Porco et 
al. 2005) supply the first evidence that impacts played an 
important role in the history of the irregular satellites of Sat- 
urn. Prior to the arrival of Cassini, Nesvorny et al. (2003) 
predicted that collisional removal between the population 
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Figure 6. Samples from Phoebe, Hyperion and lapetus, calibrated through the RC15 (2005) sensitivity function and projected on the 
different planes representing, respectively from top to bottom, the band depths of: the CH overtone signature at 1.75 fira, the 2.44 ^m 
feature, the aromatic CH stretch at 3.25 fim, and the CH stretch at 3.41 fim, with respect to the 1.51 fim water ice band depth. Left 
column: association with the satellites. Right column: homogeneous types found by the G-mode analysis; a black triangle marks the 
barycentre of the classes. 
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Figure 7. Samples from Phoebe, Hyperion and lapetus, calibrated through the RC15 (2005) sensitivity function and projected on the 
different planes representing, respectively from top to bottom, the band depths of: C02 at 4.26 /im, the CN-related feature at 4.42 fivci, 
the CN fundamental stretch at 4.53 fim, and the spectral continuum profile (normalised I/F at 3.55 /im versus normalised I/F at 1.82 
fim). Left column: association with the satellites. Right column: homogeneous types found by the G-mode analysis; a black triangle 
marks the barycentre of the classes. 
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Figure 8. Samples from Phoebe, Hyperion and lapetus, calibrated through the RC17 (2008) sensitivity function and projected on 
different planes representing, respectively from top to bottom, the band depths of: the CH overtone signature at 1.75 fj,m, the 2.44 
feature, and the aromatic CH stretch at 3.25 /im with respect to the 1.51 water ice band depth. Left column: association with the 
satellites. Right column: homogeneous types found by the G-mode analysis; a black triangle marks the barycentre of the classes. 



of irregular satellites should have caused significant modi- 
fications of the Saturnian system over the age of the Solar 
System. Phoebe, as the biggest member of the Saturnian ir- 
regular satellites, likely played a major role in this removal 
process (ibid). At the time, only 13 irregular satellites had 
been discovered around Saturn: estimating the colhsional 
evolution of this population, Nesvorny et al. (2003) reported 
that about 6 impacts (mostly due to Phoebe) were to be ex- 
pected over a time comparable to the age of the Solar Sys- 
tem. Using an updated sample counting 35 irregular satel- 



lites, Turrini, Marzari & Beust (2008) showed that, taking 
into account the present population, a number of impacts 
more than twice higher should be expected over a time of 
3.5 X 10^ years, i.e. the hypothesised lower limit to the age of 
the irregular satellites. Such results would imply an original 
population at least 33% higher than the present one (ibid). 
Moreover, the orbital structure and the radial distribution 
of the Saturnian irregular satellites suggests that Phoebe 
played a major role in clearing its near-by orbital region from 
once existing colhsional shards and smaller objects (Turrini, 
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Figure 9. Samples from Phoebe, Hyperion and lapetus, calibrated through the RC17 (2008) sensitivity function and projected on 
different planes representing, respectively from top to bottom, the band depths of: the CH stretch at 3.41 fim, C02 at 4.26 ^m, and the 
spectral continuum profile (normalised l/F at 3.55 ^m versus normalised l/F at 1.82 fim). Left column: association with the satellites. 
Right column: homogeneous types found by the G-mode analysis; a black triangle marks the barycentre of the classes. 



Marzari & Beust 2008; Turrini, Marzari & Tosi 2009). 
Finally, even if less clearly than the Jovian system, the or- 
bital structure of the Saturnian irregular satellites hints to 
the existence of possible coUisional families produced by the 
disruption of bigger parent bodies (Nesvorny et al. 2003; 
Turrini, Marzari & Beust 2008). 

Most impacts between Saturnian irregular satellites would 
involve pairs of counter-revolving, i.e. prograde vs. retro- 
grade, objects (Nesvorny et al. 2003; Turrini, Marzari & 
Beust 2008) . As shown by Nesvorny et al. (2003), the as- 



sociated impact velocities are about 3 km s~ . This implies 
that the dust-generating collisions in the Saturnian system 
would be as energetic as those in the Jovian system, which 
mainly involve members of the sole prograde or retrograde 
populations of irregular satellites and have impact velocities 
of the order of 1 — 2 km s"'^ (ibid). Also collisions between 
the Saturnian irregular satellites and interplanetary objects, 
likely at the basis of the formation of the putative coUi- 
sional families, would be characterised by the same range of 
impact velocities. As shown by Zahnle et al. (2003), col- 
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lisions of ecliptic comets on Phoebe would involve average 
impact velocities of about 3 km s~^. Since the impact ve- 
locity is mainly due to the encounter velocity between the 
host planet and the comet (ibid) and is little influenced by 
the differences in the small orbital velocities of the irregular 
satellites, we can assume the value computed by Zahnle et 
al. (2003) for Phoebe as generally valid for all Saturnian 
irregular satellites. 

In addition to these possible major impact events, it's been 
suggested by Krivov et al. (2002) that impacts of interplan- 
etary micrometeorites with small, atmosphereless bodies like 
the irregular satellites would supply a smaller yet continuous 
dust source. Such hypothesis is based on the data collected 
by Galileo space mission, which suggest that a dust produc- 
tion process is still ongoing in the system and is responsible 
for an enhancement of the circum-Jovian dust respect to 
interplanetary dust of about an order of magnitude (ibid). 
While the results of Krivov et al. (2002) refer exclusively 
to Jovian system, their implications are in principle valid 
for all giant planets with the caveat that smaller fluxes of 
micrometeorites are to be considered. As shown by Zahnle 
et al. (2003) for ecliptic comets, in fact, the impact rate 
on Phoebe is on average half than that on Himalia, the 
biggest Jovian irregular satellite. As showed by both theoret- 
ical modeling (see e.g. Cordelli &: Farinella (1997); Benz & 
Aspaugh (1999)) and observations (e.g. Meech et al. (2005); 
Ipatov & A'Hearn (2008) and references therein for ejecta 
speed from Deep Impact experiment), most of the produced 
dust particles would be ejected with velocities lower than 
those (i.e. v > 400 m s~^) needed to escape Saturn's grav- 
itational attraction (Turrini, Marzari & Tosi 2009), thus 
remaining on planetocentric orbits. 

While all these dust production processes were suggested 
by indirect evidences and comparative considerations with 
the Jovian system, the recent discovery of a ring of parti- 
cles around Saturn spanning over most of the orbital region 
of the irregular satellites (Verbiscer, Skrutskie & Hamilton 
2009) represents the proof that dust production processes 
acted and, by analogy with the Jupiter's gossamer rings, 
are still acting in the Saturnian system (ibid). Even if the 
information on this newly discovered disk is still limited, we 
argue that its formation is likely connected to the capture 
of Phoebe (Turrini, Marzari & Tosi 2009) and the origin 
of Phoebe's gap (Turrini, Marzari & Beust 2008; Turrini, 
Marzari & Tosi 2009). If collisional families are present 
between the Saturnian irregular satellites as suggested by 
Nesvorny et al. (2003) and Turrini, Marzari & Beust (2008), 
the disruptive collisions that generated them would supply 
additional dusty material and collisional shards which are 
longer-lived than those ejected in the orbital region crossed 
by Phoebe (Turrini, Marzari & Beust 2008). 

6.2 Transfer mechanism from the outer Jovian 
system 

Dust grains in orbit around Saturn would experience various 
perturbing effects, namely the gravitational perturbations of 
the Sun and the outer planets, the radiation pressure and, 
on longer timescales, the PR drag. While particles smaller 
than 1.5 jim would be rapidly ejected from the Saturnian 
system due to the effects of radiation pressure as reported 
by Verbiscer, Skrutskie & Hamilton (2009), bigger particles 



(i.e. > 3.5 ^m, ibid) would likely survive long enough for 
PR drag to act. Following the work by Soter (1974) and 
Burns et al. (1996), the orbital evolution of such dust par- 
ticles would cause them to spiral towards Saturn. If the dust 
grains where originally located outside the orbit of lapetus, 
during their radial migration they intersect the orbit of the 
satellite. During their inward drift, dust grains moving on 
retrograde orbits would impact lapetus on the leading hemi- 
sphere due to their counter-revolving motion, while the same 
argument does not apply to prograde grains. Moreover, as 
a consequence of the counter-revolving motion, retrograde 
dust particles would experience a higher frequency of close 
encounters with the satellite, suggesting that lapetus should 
be more efficient in collecting them than those in prograde 
motion. This discrepancy in the capture efhciency could jus- 
tify the asymmetry in the distribution of the dark material 
covering the surface of the satellite. Such mass transport 
mechanism is probably still active in the Solar System, as 
shown by the results of Galileo mission we mentioned above. 
The data supplied by Galileo spacecraft, moreover, indicated 
that a significant fraction of the collected dust grains were 
on retrograde planetocentric orbits (Krivov et al. 2002). 
To evaluate the efficiency of lapetus to collect the dust grains 
drifting towards Saturn due to Poynting-Robertson drag, 
we used a modified version of the algorithm developed by 
Kessler (1981). The original algorithm has been designed to 
study the collisional evolution of the Jovian irregular satel- 
lites by computing the impact probabilities between pairs 
of them. The algorithm takes into account the spatial and 
temporal density distributions of the two satellites, assuming 
that their main orbital elements (semimajor axis, eccentric- 
ity and inclination) are fixed in time and that in the times- 
pan considered their remaining orbital angles would sample 
uniformly all possible values (i.e. the timespan considered 
should be longer than their precession timescales). While 
the dust grains migrate inward, they go subject to the grav- 
itational perturbations of Jupiter and the Sun in the outer 
Saturnian system and of the regular satellites in the inner 
Saturnian system. As a consequence, the invariance of the 
main orbital elements cannot be assumed. 
To overcome this issue, instead of evaluating the impact 
probability between lapetus and real, migrating and evolv- 
ing dust grains, we used Kessler's method to evaluate the 
efficiency of lapetus in sweeping different regions of the or- 
bital phase space. We considered a radial region comprised in 
the range [1.25xl0~^ - 2.38x10"^] AU from Saturn (i.e. the 
radial region where orbits with different eccentricity values 
can intersect the one of lapetus) divided into 20 concentric 
rings. In this region we computed the probability of collision 
of lapetus with synthetic massless particles whose eccentric- 
ity and inclination^ values vary in the ranges [0-0.9] and 
[0°-60°] for prograde orbits or [120°-180°] for retrograde or- 
bits. The sampling steps assumed were respectively 0.1 and 
2°. The time over which to integrate the impact probabil- 
ity for each particle has been assumed equal to the time 
needed to cross the radial ring where it is initially located. 
This migration time is computed through the relationship 
for planetocentric Poynting-Robertson drag (see p. 430 from 



^ The inclination of the particle's orbit is measured respect to 
the planet's orbital plane about the Sun 
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De Pater & Lissauer (2001)): 

da a 5 + cos^i* 

dt tpr 6 

where tpr is a characteristic decay time given by 
1 



530^ yr, 



(2) 



(3) 



/3 is the ratio between forces due to the radiation pressure 
and the gravity of the Sun, a is the semimajor axis of the 
particle, i its inclination respect to the orbital plane of 
the planet about the Sun, tq is the Sun-Saturn distance 
and rAu is the same distance expressed in astronomical 
units. Inclination and eccentricity are not influenced by PR 
drag (ibid). As our template orbit for lapetus we used the 
mean orbital elements supplied by the JPL Solar System 
Dynamics website^. However, since the inclination referred 
to the local Laplace plane while our reference frame is the 
orbital plane of Saturn, we recomputed the mean inclination 
of the satellite (i=10.152°) from the results of Model 2 of 
Turrini, Marzari & Beust (2008). 

For our numerical setup and the grain sizes we considered 
(see section 6.3), the migration times are always longer 
(usually by at least an order of magnitude) than the 
precession timescales of the bodies involved: this allowed us 
to apply Kessler's method to each radial region. To estimate 
the impact probability of real dust grains, we computed the 
cumulative impact probability over the whole radial path 
of particles sharing the same eccentricity and inclination 
values. Globally, we considered 6510 prograde orbits and 
6510 retrograde ones: once integrated over the radial path, 
the number of orbits accounted for is 310 for each of the 
two cases in inclination. 



6.3 Dust capture efficiency 

We applied our modified model to perfectly absorbing grains 
of different sizes, which are characterised by different migra- 
tion timescales through the dependence of /3 by the particle 
radius (see e.g. p. 35 from De Pater & Lissauer (2001)). We 
considered the following values for the size of the grains: 0.1, 
1, and 10 fim. Such range of values was selected on the basis 
of comparative considerations with the results obtained by 
Krivov et al. (2002) for the Jovian system. A recent esti- 
mate performed by Verbiscer, Skrutskie & Hamilton (2009) 
for the Saturnian system suggests more strict constraints 
(i.e. d > 3.5 nm) to the size of the grains which could sur- 
vive long enough to be affected by PR drag. The size range 
we considered overlaps with the one indicated by Verbiscer, 
Skrutskie & Hamilton (2009), therefore we decided to keep 
our original one since it emphasizes the changes in the cap- 
ture efficiency as a function of the size of the grains. 
The estimated sweeping efficiency of lapetus, expressed 
through number of collisions during the crossing of the whole 
radial path, is plotted in Figs. 10 and 11 for prograde and 
retrograde dust grains respectively as a function of the ec- 
centricity and the inclination values. While lapetus has a 
capture efficiency of 1 or 100% over those regions of the 
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Figure 10. Dust capture efficiency of lapetus for increasing sizes 
of the grains (from top to bottom: 0.1, 1 and 10 microns) orbiting 
Saturn on prograde orbits. Capture efficiency is expressed as the 
number of impacts over the migration time across the whole radial 
path. Dust grains ejected by prograde satellites would likely reside 
in the region of the e — i plane (30° < i < 60° and e > 0.1) where 
lapetus' capture efficiency is lowest. 



considered phase space characterized by more than 1 colli- 
sion, we preferred to express the capture probabilities as a 
function of the number of impacts in order to emphasize the 
different statistical weights of our estimates (i.e. 3 collisions 
is statistically more significant than 1 collision). 
As can be seen from Fig. 10, the prograde dust grains cap- 
tured more efficiently are those with very low eccentricity 
values (e ~ 0): for grains bigger than 10 fim, the capture 
probability is of the order of 100% yet, once the dynamical 
features of the prograde irregular satellites are accounted 
for, this dynamical class of particles is physically unrealistic. 
There is a second island of higher efficiency at high eccen- 
tricity values (e > 0.4) for those particles having inclination 
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Figure 11. Dust capture efficiency of lapetus for increasing sizes 
of the grains (from top to bottom: 0.1, 1 and 10 microns) orbiting 
Saturn on retrograde orbits. Capture efficiency is expressed as the 
number of impacts over the migration time across the whole radial 
path. Dust grains ejected by Phoebe and the other retrograde 
satellites would reside in the region of the e — i plane (150° < i < 
180° and e > 0.1) where lapetus' capture efficiency is highest. 



values similar to that of lapetus but the capture probabil- 
ity is lowest in the region of phase space (e > 0.1 and 30° 
< i < 60°) populated by the prograde irregular satellites. 
Set aside for the high-end tail of ejection velocities, the bulk 
of the impact-generated dust grains would populate the low 
efficiency region of the e — i plane: as a consequence, the 
capture efficiency of lapetus would be of the order of a few 
0.1% for 1 /im grains and would never exceed 20-30% even 
for 10 /xm grains. 

The situation is different for retrograde grains (see Fig. 11). 
Again, the highest capture probabilities are in the low eccen- 
tricity (e ~ 0), physically unrealistic region. For retrograde 
grains, however, the island with relatively high capture prob- 



ability overlap the region of the e — i plane populated by the 
retrograde irregular satellites (i.e. e > 0.1 and 150° < i < 
180°). As a consequence, 1 /im grains are characterised by 
capture probabilities ranging between 10-40%, while for the 
majority of 10 /im grains the probabilities are of 100%. 

6.4 Mass transfer from Hyperion 

Before discussing the conclusions of this work, we would like 
to review the mass transfer scenarios linking Hyperion to the 
dark material on lapetus in light of the results of Cassini 
mission. Prior to the arrival of Cassini to Saturn, several 
studies (Matthews 1992; Marchi et al. 2002; Dobrovol- 
skis & Lissauer 2004) evaluated the transfer efficiency for 
collisionally-generated material from Hyperion to lapetus. 
Matthews (1992) assumed narrow, conical-shaped clouds of 
ejecta and estimated a single-passage transfer efficiency of 
10"^ A more detailed and realistic evaluation of the same 
scenario performed by Marchi et al. (2002) raised the trans- 
fer efficiency to about 20-40%, depending on the characteris- 
tics of the break-up event. However, the authors emphasised 
that the assumptions on the ejection direction were made to 
study the transfer efficiency of those fragments most likely 
to reach lapetus. They also noted that such particles would 
likely represent only a fraction of the ejecta cloud generated 
by catastrophic disruption events and that, if isotropic ejec- 
tion is assumed, the mean transfer efficiency lowers to 0.4%. 
Such value is in agreement with the one found by Dobrovol- 
skis & Lissauer (2004) in studying the impact probability of 
fragments ejected by Hyperion against the other Saturnian 
satellites. 

Matthews (1992) estimated the mass of dark material re- 
quired to cover the Cassini region on lapetus to a depth 
of 1 km to be about 3x10^^^ g (i.e. about half the mass of 
Hyperion, which is about 5x10^^^ g as reported by Thomas 
et al. (2007)). The assumption on the depth of the dark 
material followed from the requirement that no subsequent 
impact should excavate the whole layer of dark material and 
expose the bright material underneath. As we said in the in- 
troduction, radar measurements of the thickness of the dark 
material layer performed by Cassini constrained its depth 
between a few decimetres to about one meter. As a conse- 
quence, the amount of material needed to cover the lead- 
ing hemisphere of lapetus is at least 3 orders of magnitude 
smaller than previously thought, i.e. approximately a few 
lO'^* g or lower. 

We can estimate the cratering event needed to roughly sup- 
ply this amount of material by inverting the formula for the 
volume of a simple, bowl-shaped crater 

V' = -^(y + ^j (4) 

where h is the depth of the crater, d its diameter and V the 
volume. If we assume the depth-to-diameter ratio estimated 
for Hyperion by Thomas et al. (2007), 0.21±0.05, we have 
d ~ 5h. We then express V as Mp~^ where p=544 kg m~'^ 
(ibid). We assume for M the value computed by Matthews 
(1992) corrected for the new depth of the dark material on 
lapetus: M=3xl0^* g. By inverting eq. 4 we have 

L /24A/\3 
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This indicates that, if we assume a complete transfer be- 
tween Hyperion and lapetus, the needed amount of material 
would be supplied by excavating a crater with depth h=8.11 
km and diameter d=40.55 km, which is comprised in the 
crater size distribution observed on Hyperion (Thomas ct 
al. 2007). However, we need to take into account the effects 
of the transfer efficiency. If we assume the most favourable 
single-event scenario, i.e. a coUimated, high velocity ejec- 
tion with a 40% transfer efficiency (the maximum estimated 
by Marchi ct al. (2002)), the resulting crater would have 
h=ll km and d=55 km. If we consider a more realistic yet 
still favourable case with a 20% transfer efficiency (Marchi 
ct al. 2002), the dimensions of the crater become /i=13.87 
km and (i=69.33 km, ncaring the high-end tail of Hyperion's 
crater size distribution (Thomas et al. 2007). It should be 
noted that, in their work, Marchi et al. (2002) assume that 
the fragments are ejected in the direction that maximises 
the transfer efficiency; moreover, these author assume a uni- 
form distribution of the ejection velocities with a maximum 
value of 1.5 km s^'^, likely overestimating the contribution 
of high- velocity ejecta. Finally, as emphasised by Thomas et 
al. (2007), Hyperion is characterised by a significant poros- 
ity, of the order of 42% if the moon is mainly composed by 
water ice and likely higher if the rock fraction is significant. 
For porosity values this high, impacts can be extremely in- 
effective in excavating high-velocity ejecta since the main 
crater-forming process is compression instead of excavation 
(Holsapple et al. 2002). As reported by Holsapple et al. 
(2002), for a target porosity of 60% only 2% of the ejecta 
achieve velocities greater than 10 m s~^. This implies that 
our estimates represent a lower limit to the size of the crater 
needed to supply the dark material on lapetus and that more 
than one cratcring event is likely necessary. 
Before concluding this review of the scenarios linking Hy- 
perion to the dark material on lapetus, we would like to 
point out that the present porosity of Hyperion could be 
the by-product of a past major coUisional event and that the 
parent body of the satellite could have been characterised 
by a significantly lower porosity. Under this assumption, if 
the mass of the parent body was about 10% higher than 
Hyperion's present mass and such excess mass was collision- 
ally removed, in principle Hyperion could supply the right 
amount of material observed on lapetus. 

6.5 Comparative discussion of the mass transfer 
scenarios 

As we previously said, due to its high porosity the present- 
day Hyperion is an ineffective source for the dark material 
coating the leading hemisphere of lapetus. A most viable 
scenario connecting Hyperion to the dark material implies a 
significantly lower porosity of Hyperion's parent body and 
the collisional removal of about 1/10 of the present mass 
of the satellite. However, the excavated material would be 
delivered in the form of collisional shards and fragments, 
thus leaving open the issue of how to obtain the uniform 
dark blanket we now observe on lapetus. 
As a comparison, due to its lower porosity and higher den- 
sity (Porco ct al. 2005), Phoebe would be a more efficient 
source of ejecta. A medium sized crater like Hylas, with 
its diameter d=28 km and its depth h=4.67 km computed 
through its depth to diameter ratio 1:6 (Giese et al. 2006), 



would supply 2.43x10^* g, i.e. about the estimated amount 
of dark material. As for Hyperion, not all the material 
excavated by impacts on Phoebe will be ejected as dust: 
part of it will consist of collisional shards and fragments. 
However, the results of Turrini, Marzari & Bcust (2008) 
show that those fragments too big to be influenced by 
radiation forces would likely re-impact the satellite during 
their dynamical lifetime. The ejection-reaccretion process 
would repeat also for the material excavated by these 
secondary impacts with an ejection velocity higher than 
Phoebe's escape velocity, i.e. about 100 ms~^ (Porco et al. 
2005), with an overall enhanced production of dust. 
The divergent results we obtained for the visible and 
infrared data, however, argue against the identification 
of Phoebe as the sole source of the dark material. Due 
to the likely active collisional history of the irregular 
satellites in Saturn system, a multi-source scenario like the 
one hypothesised by Buratti, Hicks & Davies (2005) is 
more plausible and the contribution of different satellites 
could explain the observed spectral differences between 
Phoebe and lapetus. Moreover, if the conjectured nature 
of Phoebe's gap (Turrini, Marzari & Beust 2008) is real, 
Phoebe coUisionally removed once-existing prograde regular 
satellites near its orbital region and reaccreted most of 
the generated collisional shards: in such scenario, the 
contribution to the dust flux of these now-extinct satellites 
should be accounted for in interpreting the spectral features 
of the dark material. 

The existence of a disk around Saturn which spans over 
the orbital region of the irregular satellites argues in 
favour of this scenario. First, it proves that collisional, 
dust-generating processes acted and are likely still acting 
in the outer Saturnian system. Second, the lower limit to 
its mass estimated by Verbiscer, Skrutskie & Hamilton 
(2009) under the conservative assumption that the disk is 
composed by 10 /im grains indicates that even moderate 
cratcring events (i.e. producing craters of about 1 km in 
diameter) would be able to supply enough material to form 
such disk. Assuming the present disk is in a steady state, 
the accumulation rate on lapetus would be of 20 /xm Myr~^ 
(Verbiscer, Skrutskie & Hamilton 2009). Depending on the 
cratcring rate on lapetus (see e.g. the timescale estimated 
by Zahnle et al. (2003) for ecliptic comets) and the real 
size distribution of the grains composing the disk, such rate 
can possibly be enough to resupply the material excavated 
on lapetus by impacts. 



7 CONCLUSIONS 

The spectral range covered by Cassini/VIMS and the avail- 
ability of data from the close flyby of lapetus allowed us 
to perform an automatic spectral classification of the sur- 
faces of lapetus, Hyperion and Phoebe at comparably high 
spatial resolution and favourable geometry. Our aim was to 
look for spectral affinities between the dark hemisphere of 
lapetus and the two satellites that have been historically 
indicated as possible sources of the exogenous material cov- 
ering it. The classification has been performed separately 
and adopting different approaches for the visual and infrared 
portions of the spectra measured by VIMS. For the former, 
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the I/Fs measured in all the available wavelengths were used 
as variables, while several diagnostic features and the spec- 
tral continuum were used as variables for the latter, further 
separating the results obtained with the application of two 
different sensitivity functions in the calibration pipeline. 
As a general remark, striking spectral associations between 
lapetus' dark material and Phoebe or Hyperion are hardly 
found in these data. In the visible range, the G-mode anal- 
ysis confirms that Phoebe is essentially grey, with a subset 
of the spectra even showing a negative slope, while lapctus 
and Hyperion are clearly reddened. However, we found no 
association between the dark side of lapetus and Hyperion, 
since the reddening appears stronger on Hyperion than on 
lapetus while the dark material on lapetus shows a lower 
albedo than Hyperion. Moreover, any association in the vis- 
ible range would not be conclusive, as similar photometric 
correlations have been pointed out also with other small 
dark satellites moving on both retrograde and prograde or- 
bits (Buratti, Hicks & Davies 2005). In the near infrared 
range the correlation appears generally stronger between la- 
petus and Phoebe, although significant variability of some 
compounds is found on all the three satellites. We observed 
a clear correlation of most non-ice features with water ice in 
the case of Hyperion, which suggests that these compounds 
exist as trapping structures (e.g., fluid inclusions) into wa- 
ter ice, whereas the same kind of correlation - particularly 
for spectral features clearly correlated with dark material - 
is not present in the data on lapetus and Phoebe. Further- 
more, the spectra of Hyperion in the near infrared show a 
continuum profile significantly different from those of lape- 
tus and Phoebe. 

The evaluation of lapetus' sweeping efficiency in a PR drag 
based scenario yielded interesting results. The unusually 
high orbital inclination of lapetus naturally enhances its 
capture efficiency for retrograde dust grains and lower that 
for prograde grains. Statistically, a significant fraction of the 
retrograde dust particles impacts the satellite while migrat- 
ing inward for grain sizes greater than 1 /im. For grains of 
several microns in size, i.e. those that would not be removed 
by radiative forces as reported by Verbiscer, Skrutskie & 
Hamilton (2009), this fraction approaches unity. For pro- 
grade particles, on the contrary, the transfer efficiency is an 
order of magnitude lower and even for the biggest grains 
we considered it never exceeds 10 — 30%. Together with the 
hypothesised past history of the Saturnian irregular satel- 
lites, i.e. the formation of Phoebe's gap (Turrini, Marzari & 
Beust 2008; Turrini, Marzari & Tosi 2009), the capture of 
Phoebe (Turrini, Marzari & Tosi 2009) and the creation of 
coUisional families (Nesvorny et al. 2003; Turrini, Marzari 
& Tosi 2009), this result can explain the striking appear- 
ance of the satellite. The interpretation of ultraviolet data 
suggests that the delivery of the dark material is a recent or 
even still ongoing process (Hendrix & Hansen 2008). This 
can be naturally explained in the PR drag scenario through 
the effects of the impacts of micrometeorites on the irregu- 
lar satellites. Such micro-impacts, particularly those on the 
smaller retrograde satellites characterised by lower escape 
velocities, would supply a continuous source of dust as sug- 
gested by the data of Galileo mission on the Jovian system 
(Krivov et al. 2002). The existence of the newly discovered 
outer disk around Saturn (Verbiscer, Skrutskie & Hamilton 
2009) strongly supports this scenario. 



In conclusion, the results of our work argue in favour of 
a link between Phoebe and the dark material on lapetus 
due to Poynting-Robertson drag of dust particles both on 
spectroscopic and dynamical basis. While likely being the 
main actor, we argue that Phoebe was not the sole source 
which contributed in supplying the dark material. Collision- 
ally removed prograde satellites which originally populated 
the orbital region near Phoebe and the parent bodies of the 
suggested collisional families existing in the system likely 
played a role and could account for the differences observed 
in the spectroscopic data. 
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